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UNFOLDINGS OF SINGULAR HOPF BIFURCATION* 

JOHN GUCKENHEIMER, PHILIPP MEERKAMPt 

Abstract. Singular Hopf bifurcation occurs in generic families of vector-fields with two slow 
variables and one fast variable. Normal forms for this bifurcation depend upon several parameters, 
and the dynamics displayed by the normal forms is intricate. This paper analyzes a normal form 
for this bifurcation. It presents extensive diagrams of bifurcations of equilibrium points and periodic 
orbits that are close to singular Hopf bifurcation. In addition, parameters are determined where 
there is a tangency between invariant manifolds that are important in the appearance of mixed- 
mode oscillations in systems near singular Hopf bifurcation. One parameter of the normal form 
is identified as the primary bifurcation parameter, and the paper presents a catalog of bifurcation 
sequences that occur as the primary bifurcation parameter is varied. These results are applied to 
estimate the parameters for the onset of mixed-mode oscillations in a model of chemical oscillations. 
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1. Introduction. This paper presents an analysis of singular Hopf bifurcation 
i~~i in slow-fast dynamical systems with two slow variables. The setting of the problem 

^^ is as follows. Slow-fast dynamical systems are written 



ex = f{x,y) 



C where the overdot denotes a derivative with respect to a time variable t, x G i?™ are 

the fast variables, y G i?" are the slow variables, and e > is a small parameter that 
—i is the ratio of time scales. Throughout this paper we assume that there are two slow 

K*" variables and one fast variable: m = 1 and n = 2. 

The singular limit of a slow-fast system is a compound entity involving both the 
slow and fast time scales. Equation (fl]) becomes a differential algebraic equation upon 
^^ setting e = 0. This defines a vector field, the reduced system, on the two dimensional 

1^ critical manifold C given by / = at points where the projection of C onto the space 

(^ of (slow) variables y is regular. A general procedure that rescales time can be used 

^— I to extend the reduced system to the fold points of C For the desingularized reduced 

T-H system, there are two types of equilibrium points: ones that satisfy both / = and 

^ <7 = 0, and ones that lie on folds but do not satisfy g = 0. The latter are called folded 

singularities. 

The presence of multiple time scales yields new bifurcation phenomena that con- 
^ tinue to be studied extensively. Singular Hopf bifurcation in the "full" system (fTj) with 

e > and folded saddle-node bifurcation of type II [28] in the desingularized reduced 
system together comprise one of these phenomena. Folded saddle-nodes of type II 
occur at equilibrium points of the desingularized slow flow that satisfy both criteria 
for a folded singularity: they lie on folds and <? = [55] . In a one parameter family of 
systems, if an equilibrium point of the desingularized system crosses the folds of C, 
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there is a folded saddle-node of type II where it does so. In the full system, there is 
typically a Hopf bifurcation point at parameter values and state space location that 
are within 0{e) of the folded saddle-node of type II. This Hopf bifurcation is singular: 
its imaginary eigenvalues have a magnitude that is intermediate between the fast and 
slow time scales. Thus, the folded saddle-node of type II can be viewed as the singular 
limit of a singular Hopf bifurcation. 

Singular Hopf bifurcations has been observed in systems modeling action poten- 
tials in neurons as well as chemical reactions: examples are a reduced Hogdkin-Huxley 
model (see f5^ and [S^ for recent results) that has a subcritical singular Hopf bifur- 
cation, a model for chemical reactors introduced by Koper [TH] |H] [1] that has a 
supercritical singular Hopf bifurcation, and a model for an autocatalator first studied 
by Petrov, Scott and Showalter (see ^B], as well as [23], [TD] and IT for more recent 
studies) . All of these models exhibit a type of mixed mode oscillation that is directly 
related to singular Hopf bifurcation. Figures in Vo et al [3T] moreover suggest that 
singular Hopf type oscillations are also present in a model of electrical activity and 
Ca dynamics in a pituitary lactotroph, introduced by Toporikova et al [30j . 

Guckenheimer [10] presented a normal form for singular Hopf bifurcation and 
displayed some of the dynamical phenomena that occur in its unfolding. The normal 
form is 

X = {y-x^)/e 

y = z-x (1.1) 

z = — /i — ax — by — cz 

which depends upon the four parameters ii,a,b,c as well as e. New dynamical phe- 
nomena, not predicted by the Hopf theory in single time-scale systems, are found in 
this family of vector fields. This is already true in systems with one slow and one fast 
variable where there is a canard explosion [2], but the dynamics are far more com- 
plicated in systems with two slow variables and one fast variable. Guckenheimer [10] 
initiated the study of these phenomena. He observed bifurcations of equilibrium points 



and periodic orbits in the normal form ( 1.1 1, and showed that a tangency of invariant 



manifolds is present. This paper extends that work with a comprehensive analysis of 
some types of bifurcations that appear in the unfolding of the singular Hopf bifurca- 



tion normal form (1.1 ). The focus of the analysis is on phenomena local to the origin. 



i.e. those that occur at a bounded distance from the origin as e — >■ 0. In system ( 1.1 ) 



local objects or events thus approach the origin as e — t' 0, and are likely to persist 



under perturbations to system ( 1.1 1. The bifurcation diagrams and tables presented 
here are aids for the analysis of any system with two slow variables in which singular 
Hopf bifurcation occurs. This is illustrated by an example in Section 5. We describe 
our results in the framework of geometric singular perturbation theory J14j . 

The normal form ( |1.1| is invariant under a scaling transformation that eliminates 
e as a parameter: set {X,Y,Z,T) = (e~i/^a;,e~^y,e~^/^z,e~^/^i) and (A, B,C) = 
(e^/^a,e6, £-^/^c) to obtain 

X' = Y -X^ 

Y' = Z-X (1.2) 

Z' = -n -AX-BY-CZ 

Most of the analysis in this paper uses the rescaled version of the normal form. We 



present a detailed catalog of codimension 1 bifurcations occurring in system (1.2) 
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indicate where these bifurcations occur, and how these bifurcations change as addi- 
tional parameters are varied. Our strategy is to do this in three stages, first examining 
a hne in the parameter space parallel to the /x-axis, then a two dimensional slice of 
the parameter space that varies (fi, A) and finally investigating how the bifurcations 
within this slice change as {B, C) is varied. 

We regard fi as the primary bifurcation parameter and seek to determine the 



sequences of bifurcations of ( 1.2 ) that occur as fi varies while {A, B, C) remain fixed. 
This includes period-doubling bifurcations, folds of periodic orbits, Hopf bifurcations, 
homoclinic bifurcations, and a tangency of invariant manifolds, in which the two 
dimensional unstable manifold of an equilibrium in the fold region intersects the re- 
pelling slow manifold tangentially. Section 2 presents the bifurcations and a sequence 
of phase portraits for the family with {A, B, C) = (-0.05, 0.001, 0.1). 

Section 3 studies how the bifurcations identified in Section 2 depend upon the 
parameter A with {B,C) — (0.001,0.1) remaining fixed. We present a detailed two 
dimensional (/i. A) bifurcation diagram, showing a curve along which tangencies of 
invariant manifolds occur, bifurcations of the periodic orbit born in the singular Hopf 
bifurcation, bifurcations of equilibrium points occurring close to the origin, and curves 
of homoclinic bifurcations. Canard explosions occur near the homoclinic bifurcations. 
We discuss the codimension two bifurcations that occur in this restricted two param- 
eter family. Due to the numerical problems associated with the large separations of 
time-scales for very small e [5] [20], we computed the bifurcation curves of periodic 
orbits as well as the tangency curve in the bifurcation diagrams for values of e in the 
range of approximately 10~^ and lO""*. Section 4 describes how this {n,A) bifurca- 



tion diagram of system (1.2) changes as (B, C) is varied and identifies codimension 3 
bifurcations. We show some topologically inequivalent (n, A) diagrams in this section 
and give a larger portfolio in Appendix B. A table lists all the fj, bifurcation sequences 
that we found. 



Section 5 discusses a modification of the normal form (1.1) that makes the criti- 
cal manifold S-shaped. Koper [19 introduced a model of chemical oscillations that is 
equivalent to a subfamily of the modified normal form. We investigate changes that 



occur in the bifurcation diagrams of (1.1) due to the modifications of the system. 



finding that bifurcations of structures located close to the origin persist and are per- 
turbed only slightly. In the modified system, trajectories that jump from the vicinity 
of the origin approach another sheet of the critical manifold, flow along the associated 
attracting slow manifold to a fold and then jump back to the attracting slow manifold 
which comes close to the origin. This sequence of events constitutes a global return 
mechanism for trajectories to pass repeatedly near the origin. Mixed mode oscilla- 
tions 6J occur in this setting with small amplitude oscillations located in the vicinity 
of the origin and large amplitude oscillations that follow the global returns. We use 
this example to illustrate how the table of bifurcation sequences from section 4 can 
be used to estimate parameters where a tangency between the unstable manifold of 
the equilibrium and a repelling slow manifold occurs. 

Section 6 summarizes the results in the paper and discusses aspects of our anal- 
ysis that remain incomplete. We provide a table listing all bifurcation labels used in 
this paper together with the full names of the bifurcations in appendix A. Appendix 



B contains a portfolio of {fJ-,A) bifurcation diagrams for system (1.2). Appendix C 



contains supplementary data on the positions of codimension 3 bifurcations in system 



( 1.2 ). Appendix D comments on the numerical methods we used for the computation 



of bifurcation diagrams. 



2. Variations of the primary parameter. This section studies bifurcations 
and phase portraits of (1.2) that occur as /i is varied with {A, B, C) — (—0.05,0.001,0.1) 
fixed. The values of these parameters are selected so that a stable equilibrium point, 
denoted Ef, undergoes supercritical Hopf bifurcation. The distance from Ef to the 
origin is 0(e^'^). There are four bifurcations that we identify in the one parameter 
family obtained from varying fi: 

• Hopf bifurcation, 

• Tangency of the unstable manifold of an equilibrium with the repelling slow 
manifold, 

• Torus bifurcation of the periodic orbit, 

• Period doubling of the periodic orbit. 

We present phase portraits that bracket each of the bifurcations to illustrate the 
changes that occur in the dynamics of this family as fj, increases. 

Guckcnheimer [TU] gives explicit formulas for Hopf bifurcation of ( 1.2 ) and its first 



Lyapunov coefficient. Here the Hopf bifurcation occurs for /i « 0.001246. Figure 2.1 



shows the attracting and repelling slow manifolds, denoted Sa and Sr throughout the 
paper, and the one dimensional strong stable manifold (black) of the equilibrium at a 
parameter /z = close to the Hopf bifurcation. Sa is drawn in red and magenta, with 
the top (red) half tending to Ef while the bottom (magenta) half flows to X = — oo. 
The top half (cyan) of Sr consists of trajectories that flow from X = — oo, while 
the bottom half (blue) consists of trajectories that originate close to one branch of 
the strong stable manifold of the equilibrium. The existence of locally invariant slow 
manifolds of singularly perturbed systems that lie within distance 0{e) from normally 
hyperbolic critical manifolds was proved by Fenichel |9j. Attracting manifolds have 
fast foliations consisting of trajectories that approach each other in forward time; 
repelling manifolds have fast foliations consisting of trajectories that approach each 
other in backward time. As e tends to zero, the flow on the slow manifolds approaches 
the flow of the reduced system. Since we have a single fast variable, the slow manifolds 
are either attracting or repelling. Extensions of these manifolds to the vicinity of the 
fold curve intersect one another. Detailed analyses of such intersections near folded 
saddles and folded nodes were performed by Benoit [T] , Szmolyan and Wechselberger 
PS] . Wechselberger ^, Guckenheimer and Haiduc [TT], Desroches et al [3] and 
Krupa and Wechselberger [22]. Sr separates trajectories in the basin of attraction of 
Ef from trajectories that flow to inflnity in forward time. 

Figure 2.2 displays the phase portrait for /i = 0.0012715. The equilibrium point 
Ef has become a saddle-focus with a two dimensional unstable manifold W^{Ef) (red) 
whose boundary is the stable periodic orbit F (green) born at the Hopf bifurcation. 
Trajectories in the one-dimensional stable manifold ^{Ef) are plotted in black. As 
^ increases, F and ^{Ef) grow in size, getting close to the slow manifolds. At 



/i = 0.0014975, they come close to touching tangentially as illustrated in figures 2.3 



and |2.4| The multipliers of F have become complex, but they still have magnitude 



smaller than one. Note from figure 2.4 that W^{Ef) forms a "scroll" which gets much 
closer to Sr than F does. 

Sr and W^iEf) begin to intersect transversely at approximately /i = 0.00156. 
This results in bistability for VF"(-E/): some trajectories escape to X = — cx), while 
others remain in the fold region, and oscillate while approaching F, which still has 
complex attracting multipliers (see figure 2.5). 
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Fig. 2.1. Phase space when {fj,, A, B, C) = (0.0, —0.05, 0.001, 0.1). The system has an attracting 
equilibrium with a complex eigenvalue pair at {X,Y,Z) = (0,0,0). Two trajectories tending to the 
equilibrium along the strong stable eigendirections are drawn with black curves. Parts of S,- that 
approach the stable manifold of the equilibrium or escape to X = oo as t — > oo are shown in blue. 
Parts of Sa that approach the equilibrium or tend to X = —oo as t —^ oo are shown in cyan. 




Fig. 2.2. Phase space when {^i,A,B,C) = (0.0012715,-0.05,0.001,0.1). The critical manifold 
has the shape of a parabolic sheet and is plotted in gray. Subsets of Sa and Sr are shown in cyan 
and blue respectively. A black dot marks the position of Ej . Its two-dimensional unstable manifold 
W^{Ef) is plotted as a red surface, with a part close to the equilibrium removed. The singular Hopf 
periodic orbit is shown as a green curve. The stable manifold W''(Ef) of the equilibrium Ef is 
drawn with a black curve, it leaves the fold region in reverse time: one branch follows Sr and the 
other remains close to Sa for a while. 




Fig. 2.3. Phase space when {fi,A,B,C') = (0.0014975,-0.05,0.001,0.1), showing the same 
objects as figure IS. Si with the same color coding. The position of the "top" part ofW"{Ef) changed 
considerably with respect to figure \2.S\ This is to be expected, since the Sa is repelling in reverse 
time. 




Fig. 2.4. Objects in phase space crossing the plane Y = 0.5 at (fi, A, B,C) = 
(0.0014975,-0.05,0.001,0.1); a selection of trajectories in Sr are plotted using blue stars, con- 
nected with straight lines. The green dot shows the position of the stable periodic orbit's intersection 
with the plane with Y > 0. The red circles, also connected by straight lines, represent intersections 
of a single trajectory in W"(i?y) with the plane with Y > 0. 



Figure |2.6| shows a phase portrait where fi — 0.0017533. Almost all trajectories 
in W^{Ef) appear to leave the fold region, tending towards ^ = —00, with some 
trajectories following Sr for extended periods of time on their last turn before leaving 
the fold region. W'"{Ef) is no longer the boundary of the basin of attraction of F, 
which remains attracting with complex multipliers. Note the transverse intersection 
of W^"(£:/) and 5^. 

F undergoes a torus bifurcation at ^ = 0.0017829, in which its stability changes 
from complex attracting to complex repelling, and an invariant torus appears around 
it (see figure 2.7). While numerical computations suggest that at this parameter 
almost every trajectory in W'"{Ef) tends to X = —00, the position of the torus and 
the dynamics close to it may still have an impact on how many "turns" trajectories 
make before exiting the fold region. The invariant torus only exists over a short 
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Fig. 2.5. Phase space when {fi, A, B,C') = (0.0015709,-0.05,0.001,0.1), In contrast to previous 
figures, parts of the unstable manifold of the equilibrium W"{Ef) near the origin are colored in red 
and magenta, distinguishing the trajectories that escape to X = — oo from those that are attracted 
to the periodic orbit. 



parameter range: at /x = 0.0017880, trajectories in the unstable manifold of F diverge 

to X = — CXD. 

We note that the multipliers of F become real repelling for larger values of /i, 
before F undergoes a period-doubling bifurcation at approximately /i = 0.0021910 



(see figure 2.8 for a phase portrait just after the period-doubling). 



3. Bifurcation structure of a two dimensional parameter slice {B, C) = 

(0.001,0.1) . This section describes the two-dimensional bifurcation diagram of sys- 
tem (1.2), with parameters {B,C) = (0.001,0.1) fixed and (/i,^) varying. The 
diagram in figure 3.1 is typical in the sense that all other generically occurring {fi,A) 
diagrams share many of its codimension 1 and 2 phenomena and bifurcations. The 
codimension one bifurcations with the symbols that label them in figure |3.1| are 

• saddle-node bifurcation (SN) 

• (singular) Hopf bifurcation (Hopf ) 

• period-doubling bifurcations (PD) 

• fold of periodic orbit bifurcations (LPC) 

• torus bifurcations (NS) 

• tangency of invariant manifolds (T) 

• canard explosions ending in homoclinic bifurcations (S) 
The codimension two bifurcations are 

• resonances of periodic orbits (R1,R2,R3,R4) 

• degenerate homoclinic bifurcations (P) 

• zero Hopf bifurcations (ZH) 

• generalized Hopf bifurcations (GH) 
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Fig. 2.6. Phase space when {fi, A, B, C) = (0.0017533, -0.05, 0.001, 0.1). The color coding is as 
in figure \2.A but here only a part of one fundamental domain of W^{Ef) is shown. 



We briefly discuss eaeh of these bifurcations in turn. 

Saddle-node bifurcation. Recall that we assume that A and C are 0(e^/^) and 
that B is 0(e). System (1.2) can have up to two equilibria. There is only one equilib- 
rium point near the origin unless A-\- C = 0{e), in which case there are two. There 
is a saddle-node bifurcation (labeled SN in the diagrams) a,t /j, = {A + C)^/(4i3) 
and iX,Y,Z) = {Xe,X^,Xe) with X^ = {A + C)/2B. The solution curve oi ^J, = 
[A + CY /{AB) bounds a thin region of parameters where no equilibria exist. 



Hopf bifurcation and its codimension 2 degeneracies. The formulas for 
Hopf bifurcation loci in the fold region and its codimension 2 degeneracies are compli- 
cated expressions, and are best understood by disregarding terms that are of higher 
order in e. For all B ^ 0, the system has only one Hopf bifurcation that occurs 
at 0(1) distance from the origin, occurring at /i = —A^/2 — AC /2 + 0{e^/^) and 
{X,Y,Z) « {A/2, A^/ A, A/2) + 0{e^l'^), making it a singular Hopf bifurcation. The 
Hopf bifurcation is degenerate at a zero Hopf bifurcation (marked ZH in the dia- 
grams), where the Hopf equilibrium point has a zero eigenvalue in addition to the 
purely imaginary eigenvalue pair, and at two generalized Hopf bifurcations (GH), 
where the first Lyapunov coefficient vanishes. The Hopf center manifold is attracting 
for values of A greater than that of the zero Hopf bifurcation, and repelling for val- 
ues of A smaller than that of the zero Hopf bifurcation. The criticality of the Hopf 
bifurcation changes at each of these these codimension 2 degeneracies. Supercritical 
[subcritical] Hopf bifurcations at values of A greater than that of the ZH bifurcation 
give rise to periodic orbits when the Hopf curve is crossed with /i increasing [decreas- 
ing] , while at smaller values of A supercritical [subcritical] Hopf bifurcations give rise 
to periodic orbits as the Hopf curve is crossed with /i decreasing [increasing] . 
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Fig. 2.7. Phase space when {n,A,B,C) = (0.0017846,-0.05,0.001,0.1). A segment of a tra- 
jectory on the invariant torus is drawn with a dark-green line. Otherwise, the color coding is as in 
the previous figure. The periodic orbit is largely obscured by the torus. 




Fig. 2.8. Phase space when {ii,A,B,C) = (0.0022000,-0.05,0.001,0.1). A period- doubled 
periodic orbit is plotted with an olive thick line. The color coding for the other periodic orbit, the 
slow manifolds, critical manifold, and invariant manifolds of the equilibrium near the origin is as 
in the earlier figures. A part of the unstable manifold of the equilibrium near the origin not plotted 
here comes very close to the period- doubled periodic orbit. 




Fig. 3.1. {p,A) bifurcation diagram for {B,C) = (0.001,0.1). 



Period-doubling bifurcations. Numerical continuation with MatCont [7j and 
AUTO |8] show that F undergoes a period-doubling bifurcation (labeled PD in bi- 
furcation diagrams) as /i is increased. In the rescaled system (1.2 1, these bifurcations 
typically occur on periodic orbits that are 0{1) distant from the origin. Near the 
endpoints of the period-doubling bifurcation curve, where the PD curve meeds a ho- 
moclinic bifurcation, the amplitude becomes even larger and the bifurcating orbits of 
system (1.1 1 no longer shrink to the origin as e — )■ 0. This degeneracy is described 
further in the paragraph on "Degenerate homoclinic bifurcations" below. 



Fold of periodic orbits bifurcations. A fold of periodic orbits bifurcation 
curve (labeled LPC) emanates from each of the generalized (singular) Hopf bifurca- 
tions (GH). This bifurcation, in which two periodic orbits annihilate each other as 



/i decreases, is initially local to the origin in system (1.1), i.e. remains at bounded 
distance from the origin as e — >■ 0, but later becomes non-local in a canard explosion, 
and ends where the LPC curve meets a homoclinic bifurcation (see the paragraph on 
"Degenerate homoclinic bifurcations" for more details). 



Torus bifurcations. A locus of torus bifurcations (NS) emanates from the zero 
Hopf bifurcation (ZH). Associated with this are attracting invariant tori that sur- 
round the bifurcating singular Hopf periodic orbit. These tori apparently exist only 
in a thin strip of parameters to the right of the torus bifurcation curve. 



Tangency of invariant manifolds. In figure 3.1 Ef has a two dimensional 



unstable manifold in the parameter region bounded below by the saddle-node curve, 
and on the left by the Hopf curve. In part of this region, VF"(i?/) intersects Sr- This 
intersection is tangential at the tangency curve (T), where the manifolds begin to 
intersect as /i increases. Note that there are values of A where the manifolds intersect 
immediately after the Hopf bifurcation. See appendix D for details on the computa- 
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tion of the tangency curve. 

Canard explosions and homoclinic bifurcations. The periodic orbit F grows 
rapidly in canard explosions. The canard explosions consist of periodic orbits that 
have long segments lying near Sr- The small parameter ranges in which they occur are 
determined by the geometry of the slow flow as illustrated in figure |3.2[ The slow flow 



is computed by differentiating the equation of the critical manifold in system (l.ll 
with respect to t, and using the result and the equations for the derivatives of the slow 
variables to obtain equations for x and z in x and z only. After rescaling the vector 
field so that the direction reverses on the repelling sheet of the critical manifold, we 
obtain the following desingularized slow flow equations: 

, (3.1) 

z = —2x {fi + ax + bx" + cz) 

Trajectories that follow S^ can jump at any location along the manifold. The singular 



limit of trajectories of system (1.1) approach concatenations of trajectories of the 
slow flow together with segments parallel to the x-axis that begin on the repelling 
sheet of the critical manifold and end on the attracting sheet of the critical manifold. 
Since the critical manifold is a parabolic cylinder symmetric with respect to reflection 
in the (j/, z) plane, the singular limits of some periodic orbits consist of trajectories 
of the slow flow that connect points (x, z) to points {-~x,z) concatenated with the 
horizontal segment from (— a;, z) to {x,z). We call these periodic orbits with a single 
fast segment simple. We refer to their singular limits as singular cycles. 

We expect the simple periodic orbits to grow rapidly where the slow flow equations 
satisfy 

z{x, z) z{--x, z) 

x{x,z) x{—x,z) 

With the simplifying assumption fi — 0, substitution of the slow flow equations into 
the ansatz gives that when fi — 0, canard explosions are likely to occur at parameters 
satisfying 

a^ +ac + b = 0, (3.2) 



and that the trajectory segments on the slow manifolds of system (1.2 1 can approx- 
imately be parametrized by {X, Y, Z) ^ {X, X'^, AX"^). Note that e(a^ + ac + h) = 



A + AC -\- B iw the rescaling from system (l.ll to system (1.2). Numerical calcu 



lations with MatCont and AUTO show that these estimates are very accurate, and 



that the location of periodic orbits is insensitive to the small value of p.. In figure 3.1 
canard explosions occur along lines that are nearly horizontal, and end in homoclinic 
bifurcations with eq uilib ria that are far from the origin. Dashed green lines labeled S 



are drawn in figure 3.1 aX A= ^JC^ — 4B as approximations to the region where the 



slow flow has singular cycles. 

Resonances of periodic orbits. There are 1:1, 1:4, 1:3 and 1:2 resonances on 
the torus bifurcation curve (labeled Rl, R4, R3 and R2 respectively), where F has 
multipliers on the unit circle with arguments 0, 7r/4, 7r/3 and tt/2 respectively. At the 
1:1 resonance, the torus bifurcation curve ends at a fold of periodic orbits bifurca- 
tion locus in a bifurcation that is analogous to the Bogdanov-Takens bifurcation for 
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Fig. 3.2. Phase portraits for parameters {ii,A,B,C) = (2.5008 X 10"*', -0.090281, 0.001, 0.1) 
(top) and {fj.,A,B,C) = (2.5008 x 10"'', -0.088758, 0.001, 0.1) (bottom), showing the projection of 
a periodic orbit (magenta, computed in the full system with MatCont) onto the slow flow. The 
nuUclines of the slow flow as well as the fold curve are drawn using red dashed lines. Segments 
of trajectories of the slow flow are plotted in blue and cyan, the two middle ones ending at the 
origin. A reflection of the middle cyan trajectory about the Z axis is plotted as a dashed cyan line. 
The position of the equilibrium far from the origin on the repelling sheet of the critical manifold is 
marked with a dot. 



2-dimensional vector fields [5^1 [H]. At the 1:2 resonance, the torus bifurcation curve 
ends at the period-doubling curve in another bifurcation analogous to the Bogdanov- 
Takens bifurcation with symmetry |291 I15j . 



Degenerate homoclinic bifurcations. The homoclinic bifurcation is degen- 
erate at two points, labeled P in the (/i. A) bifurcation diagrams. At each of these 
points, a fold of periodic orbits curve (LPC) and a period-doubling curve (PD) meet 
the curve of homoclinic bifurcations. We believe that these codimension 2 phenomena 
are due to homoclinic orbit flips (top point) and homoclinic inclination flips (lower 
point). Unfoldings of these bifurcations were developed by Sandstede [37] (orbit flip) 
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and Homburg et al [T3j (inclination flip) and Kisaka et al. [17], [T8| (inclination flip). 
We did not investigate the points P in detail. 



4. Catalog of bifurcation diagrams in one and two dimensional param- 
eter slices. This section first describes curves of codinicnsion 1 bifurcations in (/i, A) 
bifurcation diagrams for different fixed values of {B,C). This includes information 
about the codimension 2 bifurcations at which codimension 1 bifurcation curves meet 



or end. We then complete our study of the bifurcations of system (1.2) by investi- 
gating how the {fi, A) bifurcation diagrams change as B and C vary. These changes 
occur when we encounter points of codimension larger than two and when we en- 
counter intersections of the loci of codimension one and two bifurcations that are not 
transverse. Most of the changes we find are of the latter type. Indeed, we identify 
here only a single degeneracy of the first type, which occurs when B = 0. 



Equilibrium point bifurcations. Unless B — 0, system (1.2) has a saddle- 
node bifurcation at fj, = {A + C)^/(4i3), whose properties are described in section 3. 
Similarly, there is only one curve of local Hopf bifurcations for B ^ 0, occurring at 
/I « -A'^/2 - ylC/2 + 0(e^/^). It is degenerate at a zero Hopf (ZH) bifurcation at 
A — C{B — 1). A computation of the normal form coefRcients using notation from 
Kuznetsov [53] shows that if B < 0, then s — 1, 6 < 0, and the system has homoclinic 
orbits, invariant tori, and invariant spheres, which are local in the vicinity of the zero 
Hopf bifurcation. If i? > 0, then s = 1, 9 > 0, and the dynamics in the vicinity of the 
zero Hopf bifurcation is simpler. 

The Hopf bifurcation is also degenerate at generalized Hopf bifurcations (GH), 
where the first Lyapunov coefficient vanishes. Simplifying the expression for the first 
Lyapunov coefficient by eliminating terms that are of higher order in e, we observe 
that there are generalized Hopf bifurcations approximately when A^ + AC + 2B = 
has a real solution, i.e. whenever C^ — 8B > 0. We note that equilibrium points at 
0{l/e) distance from the origin can undergo Hopf bifurcation as well. 

Periodic orbit bifurcations. Period-doubling bifurcations (PD) exist in those 
{B, C) regions where the slow flow has singular cycles. The bifurcation loci begin 
and end at points labeled P, where the periodic orbit that is to double undergoes a 
canard explosion. At these points, the period-doubling bifurcation curve meets a fold 
of periodic orbits curve, in what seems to be a tangential way. The period-doubling 
occurs as fi increases through the bifurcation parameter value, and the first period- 
doubling bifurcation may be followed by further period-doublings. Fold of periodic 
orbits (LPC) curves exist in those {B, C) regions where the slow flow has singular 
cycles. There is one fold of periodic orbits curve if no generalized Hopf bifurcation 
is present, and the bifurcation curves begin and end at meeting points labeled P, 
where they connect with period-doubling bifurcations. If two generalized Hopf bifur- 
cations are present, then the fold of periodic orbits curves begin at generalized Hopf 
bifurcations, and end at points marked P. Given {B,C), the sign of B determines 
whether the LPC curve hes to the left {B > 0) or to the right {B < 0) of the Hopf 
bifurcation curve. Canard explosions of the slow flow occur for parameters near those 
that satisfy A^ + AC + B = where the slow flow has singular cycles. These canard 
explosions typically end in homoclinic bifurcations with equilibrium points that are 
at 0(e) distance from the origin. Being non-local in many cases, we did not investi- 
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gate honioclinic bifurcations thoroughly and plot homoclinic bifurcation curves only 



in figure B.l which gives the {fJ.,A) bifurcation diagram for {B,C) = (—0.01,0.1). 

Invariant tori exist close to the torus bifurcation curves marked NS. If -B < 0, 
normal form theory predicts that these tori are unstable close to the zero Hopf bi- 
furcation [5S], and exist for values of A greater than that of the torus bifurcation 
curve. A more detailed analysis of the stability and persistence of the torus bifur- 
cations was not performed. When B < 0, the torus bifurcation curve starts at the 
zero Hopf bifurcation, and ends at a 1:2 resonance (labeled R2), where it meets the 
period-doubling curve. The torus bifurcation curve intersects no other bifurcation 
curves included in the diagrams. The defining equation that characterizes torus bi- 
furcations is also satisfied by "neutral" periodic orbits with two multipliers whose 
product is one. Consider first parameters {B,C),B > for which there is no gen- 
eralized Hopf bifurcation. A curve i(_B,c) of neutral periodic orbits emanates from 
the Hopf bifurcation. This curve is shown in figure |B.10[ but is omitted in all other 
bifurcation diagrams. Close to the zero Hopf bifurcation, the multipliers are real, and 
the periodic orbits are saddles. The curve i(B,c) can intersect the period-doubling 
bifurcation curve in 1:2 resonances or meet a fold of periodic orbits bifurcation curve 
tangentially in 1:1 resonances where the multipliers become complex. Both events can 
occur more than once in the same (fi, A) diagram, due to bending of the curve i(B,c)) 
see e.g. figure |B.8| for two 1:2 resonances, or figure |B.10[ If there is a generalized 
Hopf bifurcation and i? > 0, then there is typically one curve of torus bifurcations, 
beginning at a fold of periodic orbits bifurcation at a 1:1 resonance, and ending at an 
intersection with a period-doubling bifurcation curve at a 1:2 resonance. Independent 
of whether a generalized Hopf bifurcation is present or not, when B > and C > 0, 
as the absolute value of C decreases towards the parameter C = C{B) where an isola 
vanishes, the curve Li^.c) tends to move along the families of periodic orbits, from 
the singular Hopf bifurcation towards and beyond the period-doubling bifurcation. 
See figures |B.3| to |B.10| for a sequence of diagrams illustrating this trend. If i? > 
and C < 0, then the families of periodic orbits are traversed away from the Hopf 
bifurcation as C increases towards the point where the isola vanishes. 

Tangency of invariant manifolds. We now discuss intersections of W'^{Ef) 
with Sr- Regions in parameter space where these manifolds intersect are separated 
from those where they do not intersect by tangencies of invariant manifolds as well as 
by three codimension 1 bifurcations: saddle-node bifurcations, Hopf bifurcations, and 
fold of periodic orbits bifurcations. The tangency curve is in many cases very close 
to torus bifurcations or period-doubling bifurcations, to the point of being indistin- 
guishable from these bifurcation curves in some of the (/z, A) bifurcation diagrams. 
Nonetheless, the tangency of invariant manifolds is a bifurcation in its own right, and 



does not coincide with period-doubling or torus bifurcations in system (1.2) except at 
isolated points of the (/j,, A) diagrams. 

The tangency curve is situated in the (//, A) diagrams in different ways for _B > 
and B < Q. li B < 0, Ef leaves an 0(1) neighborhood of the origin in system 



(1.2) when A + C ~ 0(e) with A + C < 0. The tangency curve begins at the 
saddle-node curve at a value of A slightly more negative than that of the zero Hopf 
bifurcation, indistinguishably close in most of the diagrams of this paper. In this 
region of parameter space, W"{Ef) may approach the one-dimensional stable manifold 
of the second equilibrium in the vicinity of the origin and then leave the fold region, or 
can "get caught" in invariant structures around the two equilibrium points involved in 
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the zero Hopf bifurcation, cf. [23] figure 8.22. The tangency curve may intersect the 
the PD curve, but systematic numerical computations suggest that the tangency curve 
does not lie visibly to the right of it at the resolution of the (/i. A) bifurcation diagrams 
included in this paper. The tangency curve ends at the top point P, where the 
boundary of the region with an intersection of invariant manifolds begins to coincide 
with the fold of periodic orbits bifurcation curve, and the intersection begins to be 
transverse at the boundary parameters. 

If -B > 0; there is a saddle-focus in the vicinity of the origin in the region R(b.c) 
that lies to the right of the Hopf curve and above the saddle-node curve. If now 
C < 0, then W'^{Ef) and Sr intersect transversely in the entire region R(^b,c)- If 
C > 0, then in some parts of the region Rib.c) the two invariant manifolds do not 
intersect. More precisely, there is a tangency curve, beginning on the Hopf bifurcation 
curve, very close to the torus bifurcation curve or at a larger value of A, lying to the 
right of the Hopf curve until it meets the Hopf curve again at a larger value of A, 
very close to the period-doubling curve or at a smaller value of A. The tangency 
curve can intersect the torus and period-doubling curves for large enough values of 
C = C{B), and there may be interactions between W'^{Ef) and the attracting tori for 
(/I, A) slightly below the torus bifurcation curve. Systematic numerical computations 
suggest that the tangency curve does not lie visibly to the right of the period-doubling 
curve at the resolution of the (/x. A) bifurcation diagrams included in this paper. It 
never lies much to the right of the torus bifurcation curve either, see figure |B.6| for a 
case in which the tangency curve lies visibly to the right of the torus bifurcation curve. 

Remark: local and non-local bifurcations. The bifurcations described above 
are, with two exceptions, local to an 0{e^^^) neighborhood of the origin when scaled 
back to the system ( |1.1[ ). The first exception is the saddle-node bifurcation (labeled 
SN), which is local in the vicinity of the zero Hopf bifurcation (ZH), but is generally 
non-local. Note also that the range of values of A for which saddle-node bifurcations 
occur for {B, C) fixed and ijl = 0{e) is 0(e), while most other codimension one bifur- 
cations are of interest over an 0(e^/^) range of A. The second exception is that the 
periodic orbit V can become non-local at parameters close to the horizontal dashed 
green lines in the (/i, A) bifurcation diagrams where the slow flow has singular cycles. 
In particular, period-doubling (PD) and fold of periodic orbits (LPC) bifurcations 
become non-local at the points marked P, and torus bifurcations (labeled NS) where 
a torus forms around a periodic orbit can become non-local if they are in the vicinity 
of the green dashed lines. 



Global bifurcation structure. System (1.2) is highly degenerate when B = Q. 
It has only one equilibrium point, which is in the fold region unless A + C = 0{e). 
Asymptotic formulas for the first Lyapunov coefficient of the singular Hopf bifurca- 
tion when B y^ suggest that the generalized Hopf bifurcation passes through the 



zero Hopf bifurcation at i? = in a non-generic manner, see figures B.l B.2 and 



B.4 3.1 B.3 for sequences of {fi, A) diagrams in which B becomes small relative to C. 
Also, as the generalized Hopf bifurcation passes through the zero Hopf bifurcation, 
the normal form type of the zero Hopf bifurcation changes [53j . We conjecture that 



the addition of higher order terms to the system ( 1.1 ) can produce qualitative changes 



in its bifurcation diagram, but do not pursue that issue here. 

We next list phenomena that produce changes in the (/i. A) bifurcation diagrams 
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due to non-transverse intersections of these slices with bifurcation manifolds of codi- 
mensions one and two: 

• folds of generalized Hopf bifurcations 

• folds of the curves where the slow-flow has singular cycles 

• tangential intersections of the period-doubling bifurcation locus and the Hopf 
bifurcation loci 

• tangential intersection of the torus bifurcation locus and the Hopf bifurcation 
locus 

• resonant zero Hopf and generalized Hopf bifurcations 

• endpoint of the tangency curve crosses the torus bifurcation locus 

• endpoint of the tangency curve crosses the period-doubling bifurcation locus 

• the tangency curve intersects the period-doubling bifurcation locus tangen- 
tially 

• the tangency curve "attaches" to the Hopf curve. 

Figure 4.1 displays 16 regions in the {B,C) parameter space of system (1.2). 
These regions constitute a partitioning into regions corresponding to classes of topo- 
logically equivalent {fJ-,A) bifurcation diagrams, with two caveats: first, regions la 
and lb as well as VHIa and VHIb have topologically equivalent (/i. A) diagrams, we 
merely distinguish them them for symmetry reasons. Changing the sign of parameter 
C corresponds to reflecting bifurcation curves in (/i. A) diagrams about the /Lt-axis, 
and to changing the time orientation. Under this transformation, the tangency curve 
transforms into a curve corresponding to a tangential intersection of ^{Ef) with 
Sa- The tangency involving Sr is present in quadrants I, H and HI of figure 4.1 
whereas the tangency involving Sa is present in quadrants H,IH and IV of figure 4.1 



Appendix B contains at least one representative (/i. A) bifurcation diagram for each of 
the 8 regions of {B, C) parameter space in which C > 0. Secondly, certain codimen- 
sion three events related to the torus bifurcation and its degeneracies are not included 



in figure 4.1 



The degeneracy at i? = (thick solid line in figure 4.1) separates regions in 
{B, C) parameter space with vastly different (/i. A) bifurcation diagrams. For C fixed, 
ifijA) diagrams become degenerate as |i3| tends to and the zero Hopf bifurcation 
approaches generalized Hopf bifurcations. There is only one topological type of (/i. A) 
diagram for B < 0. The situation is more complicated for B > 0. What follows is a 
description of the topological changes in the (/i, A) diagrams occurring at the curves 
separating the regions where B > 0. We only treat the first quadrant of figure 4.1 a 



description for the fourth quadrant may be obtained using the symmetry C h^ — C. 



Note that no curves in figure 4.1 intersect, except maybe for < i? <C 10 . 

Regions Ha and IHa are separated by curve E, which consists of {B, C) values 
where the surface of generalized Hopf (GH) bifurcations is tangent to the (/i, A) slice 
of the parameter space. The formula for the first Lyapunov coefficient indicates that 
the curve E is tangent to the parabola 8i3 — C^ at the origin of the [B, C) plane. On 
this curve, there is a single generalized Hopf (GH) bifurcation of the (/i. A) bifurcation 
diagrams. To the left of this curve, there are two generalized Hopf bifurcations; to 
the right there are none. Curves F, G, H and I all correspond to events relating to 
the tangency curve or to bifurcations of periodic orbits. These curves were found by 
inspecting appropriate sequences of (/i. A) diagrams. A superset of the data used can 
be found in appendix C. Regions IHa and IVa are separated by curve F, where the 
endpoint of the tangency curve crosses the torus bifurcation locus. In region IHa, the 
tangency curve and torus bifurcation curve intersect, while they are disjoint in region 
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Fig. 4.1. Regions in {B,C) parameter space. Solid and dashed lines indicate the positions of 
codimension 3 phenomena that separate regions in {B, C) parameter space yielding different corre- 
sponding {fi,A) bifurcation diagrams. 



IVa. The tangency curve coalesces with the Hopf bifurcation curve at the boundary of 
regions IVa and Va, represented by curve G. Regions Va and Via are separated by the 
curve H, which consists of (B, C) values where the torus bifurcation locus intersects 
the Hopf bifurcation locus tangentially: to the right of H, the torus bifurcation locus 
and the Hopf curve no longer intersect. In region Via, the period-doubling bifurcation 
locus intersects the Hopf bifurcation in two points. As B is increased, these two points 
coalesce when the two bifurcation loci intersect tangentially on the curve I, after which 
they no longer meet in regions Vila and Villa. Regions Vila and Villa are separated 
by the curve J where the two (/i. A) curves of homoclinic bifurcations coalesce as B 
increases. The curve J is approximated by the quadratic equation C^ — AB = 0. 

Torus bifurcation and the global bifurcation structure. For B < 0, tak- 
ing account of the torus bifurcation curve and its degeneracies (resonances) does not 
yield more topologically different (/i, A) bifurcation diagrams. For B > 0, the division 
of {B,C) parameter space into regions with topologically equivalent (/i,A) diagrams 
gets significantly more complicated and is not presented here. The paragraphs of sec- 
tion 4 headed "Periodic orbit bifurcations" give a general description of the positions 
of torus bifurcation curves relative to other bifurcation curves. Data on the position 
of the crossing of a 1:2 resonance with the Hopf curve can be found in appendix C. 



One dimensional bifurcation diagrams. Table 4.1 lists the sequences of lo- 
cal bifurcations shown in the (/i. A) diagrams of this paper that occur as the primary 
bifurcation parameter /i is increased. We list only sequences occurring well within 
the {B, C) regions we have identified. There could be some additional sequences, e.g. 
extremely close to the zero Hopf bifurcation, or where \B\ << C^, i.e. where the zero 
Hopf, generalized Hopf and P bifurcations occur at very similar values of A. Two 
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Sequence # 


Sequence description 




[B, C) region present 


1 


Hsup 




la, Ha, Ilia, IVa, Va, lb, lib, Illb, IVb, Vb 


2 


Heup - (SN) 




Ha, Illa, IVa, Va, lib, Illb, IVb,Vb 


3 


(SN) - Hsup - T ± PD 




lb 


4 


(SN) - H,up - T ± NS - 


PD 


la 


5 


H,up - LPC 




la, lb 


6 


Heup - T ± PD 




la. Ha, lb 


7 


Hsup - T ± NS - PD 




Ha 


8 


LPC - Hsup - NS - PD 




lib, HIb 


9 


LPC - Hsup - PD 




lib, HIb 


10 


LPC - Hsup - T ± PD 




Ila 


11 


LPC - NS - Hsup - PD 




lib, HIb 


12 


LPC - PD - Hsup 




lib, HIb, IVb 


13 


Hsub 




la. Ha, Ilia, IVa, Va, lb, lib, Illb, IVb, Vb 


14 


Hsub - (SN) 




Ha, Illa, IVa, Va, lib, HIb, IVb, Vb 


15 


(SN) - Hsub - PD 




la 


16 


(SN) - Hsub - NS - PD 




lb 


17 


Hsub - LPC 




la, lb 


18 


Hsub - PD 




la, lb, lib 


19 


Hsub - NS - PD 




lib 


20 


LPC - Hsub - NS - PD 




Ha, Hia 


21 


LPC - Hsub - T ± NS - 


PD 


Ha, Ilia 


22 


LPC - Hsub - T ± PD 




Ila, Ilia 


23 


LPC - Hsub - PD 




Ila, HIa 


24 


LPC - NS - Hsub - T ± PD 


Ha, HIa 


25 


LPC - PD - Hsub 




Ha, HIa, IVa 



Table 4.1 

List of bifurcation sequences occurring in system ( |1.2| l as the main bifurcation parameter fj, 
is increased. Each sequence is described in terms of the codimension 1 bifurcations involved, fx 
increasing. The sign "-" separates codimension 1 bifurcations occurring at different values of /i. 
A "± " sign appears between two bifurcations where the first bifurcation listed either occurs before 
the second or indistinguishably close to the second, at the resolution of the diagrams in this paper. 
Saddle-node bifurcations are listed in parentheses, since they are only local close to the zero Hopf 
point. 



branches of periodic orbits emanate from LPC bifurcation points. If a sequence begins 
with an LPC bifurcation and also contains a torus bifurcation or period-doubhng, then 
one branch of periodic orbits emanating from the LPC bifurcation has a torus bifur- 
cation and period-doubhng, while the other one ends in a singular Hopf bifurcation. 
W^lEf) and Sr intersect after the tangency in the bifurcation sequences. 

5. An Example: the Koper Model. This section demonstrates how our nor- 
mal form analysis can be used in investigations of systems undergoing singular Hopf 
bifurcation. We use a vector field that was studied first by Koper T5] as a model of 
mixed mode oscillations in chemical systems: 

eix = ky — x^ + 3x— X 

y ^ X — 2y + z (5-1) 

Further discussion of this model can be found in [^, [TB] as well as Desroches et 
al. [5] and references cited there. Using affine coordinate changes and a time rescaling 
(see [6]), the Koper vector field can be written in the following form: 

X — {y — x^ ~ x^)/e 

y ~ z — X (5-2) 

i = — /i — ax — by ~ cz 
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Note that 



the system (5.2 ) has one more parameter than the system (5.1 ), so the Koper 



vector field is a subfamily of (5.2), and 



the system (5.2 1 is a variant of the singular Hopf normal form ( 1.1 ) in which 



the parabolic critical manifold has been replaced by a cubic critical manifold 
[6] by the addition of a single, higher order term to the equation. 



With the scaling of parameters and coordinates used to obtain system (1.2) from 
system (|1.1[), (5.2) transforms to 



X = {Y -X^ -e^/^X^) 

Y ^ Z -X 

Z = -1^1- AX - BY - CZ 



(5.3) 



It is apparent that the system (5.3) is a small perturbation of the system (5.2) in 



a bounded region of (X,Y,Z) space. Figures 5.1 and B.l show (/i, A) bifurcation 
diagrams for systems (5.3) and (|1.2[) at the same values of {B, C) = (- 



-0.01,0.1). The 

similarity of these two figures supports the use of the normal form ( |1.1[ ) to study local 
bifurcations near a singular Hopf bifurcation. Note that the cusp of the saddle-node 
curve in figure 5.1 is not local: it tends to infinity as £ — > 0. 




Fig. 5.1. Bifurcation diagram for cubic singular Hopf normal form, B = —0.01, C = 0.1. 

Consider the Koper model with £2 = 1- When additionally £i =0.1 and k = —10, 
there is a supercritical singular Hopf bifurcation at approximately A « —7.670. Con- 
tinuing this family of periodic orbits with varying lambda using AUTO, we find that 
the emerging periodic orbit undergoes a period-doubling bifurcation at A « —7.461, 
and the "undoubled" periodic orbit undergoes a fold of periodic orbits (LPC) bifur- 
cation at A « —6.235. The periods of the periodic orbits at the bifurcation points 
are approximately 0.64, 0.82 and 1.67 respectively. If r denotes the initial period of 
the periodic orbits emerging from the singular Hopf bifurcation, local bifurcations of 
periodic orbits typically occur with periods less than 1.5r. This suggests that the 
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LPC bifurcation is not local. Figure |5.2| shows that indeed, the shape of the periodic 
orbit at the LPC bifurcation follows the cubic shape of the critical manifold, rather 
than staying close to the fold curve at which the singular Hopf bifurcation occurred. 




1.8 
1.6- 
1.4- 
1.2- 



0.6 - singular Hopf 
(supercritical) 

0.4- 
0.2 

o' ' 



Fig. 5.2. Periodic orbits emanating from a supercritical Hopf bifurcation for system dSTTf 
and ei = 0.1, fe = —10. Observe that the periodic orbits follow the cubic shape of the (parameter- 
dependent) critical manifold more as A increases. 



Among the sequences with a supercritical Hopf bifurcation, only 3, 4 and 6 to 12 
also have a period-doubling. Note that sequences 8 to 12 have a period-doubling bifur- 
cation and a supercritical Hopf bifurcation on distinct branches of periodic orbits, and 
that sequences 4 and 7 contain torus bifurcations. Thus, only sequences 3 and 6 match 
the continuation data. This suggests that there is a tangency of invariant manifolds 
either between the supercritical Hopf bifurcation and the period-doubling bifurcation 
or very close to the period-doubling bifurcation: indeed, numerical calculations show 
that the tangency occurs at A ~ —7.539. 

Due to the S-shape of the critical manifold in the Koper model, trajectories in the 
unstable manifold in the vicinity of the fold region that escape from the fold region 
can return to the fold region. As a result, the tangency of invariant manifolds can 
mark the onset of mixed-mode oscillations: a trajectory approaches the equilibrium 
point near the fold region, then follows the equilibrium's two-dimensional unstable 
manifold, leaves the fold region, and approaches the equilibrium once more along its 
stable manifold. The Koper system has mixed mode oscillations at A = —7.50, il- 



lustrated in figure 5.3 Thus local analysis of system (1.2 1 led to information about 



the Koper model as an alternative to applying numerical continuation tools directly 
on the Koper model. This approach can be used to explore the effects of varying 
additional parameters in the Koper model as well as to study other systems in which 
singular Hopf bifurcation occurs. 
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Fig. 5.3. Mixed mode oscillations in system (jHTTJ with ei = 0.1, £2 = 1, fe = —10, A = —7.50; 
the left panel shows the time-series of the x-coordinate of an MMO trajectory, the right panel shows 
the same trajectory in the xy-plane. Ef is marked with a black dot, the critical manifold is drawn 
with a dashed black line. 



6. Discussion. This paper includes a comprehensive analysis of of bifurcations 
that occur in the normal form (1.1) for a singular Hopf bifurcation. In particular, 



we analyze bifurcations of equilibrium points, periodic orbits originating from these 
equilibria and "simple" tangencies between the unstable manifold W^^Ef) of an equi- 
librium Ef and a repelling slow manifold Sr of the system. We approach the study 
of these bifurcations in a hierarchical fashion, designating the parameter fi as the 
primary bifurcation parameter and then performing analyses of parallel two dimen- 
sional slices throughout the four dimensional parameter space of the system. While 
our analysis is comprehensive, it is hardly complete. We point to four aspects of the 
bifurcations of this system for further study. 



In some parameter ranges, the system (1.1 1 has complicated invariant sets in 



0{y/e) neighborhoods of the origin. Torus and period-doubling bifurcations of periodic 
orbits give rise to invariant sets that undergo further bifurcations. We have not 
pursued the study of these invariant sets or their bifurcations. Trajectories in these 
sets typically make multiple oscillations when projected to the (x, y) coordinate plane. 
Similarly, there are more complex intersections of invariant manifolds than those we 
have investigated in this paper. The three two dimensional invariant manifolds that we 
have studied are the attracting and repelling slow manifolds and the unstable manifold 
of an equilibrium point. In some parameter regimes, the attracting and repelling slow 
manifolds have multiple intersections that constitute maximal canards of a folded 
node point of the system [H [11] [32] ■ These manifolds and their intersections bound 
rotational sectors of trajectories that make different numbers of oscillations as they 
pass through the folded node region. The unstable manifold of the equilibrium can 
also become entangled with the repelling slow manifold and have multiple intersections 
that separate regions of the unstable manifold that flow to x = — oo after making 
different numbers of oscillations. In systems like the Koper model where trajectories 
jumping from the vicinity of a singular Hopf bifurcation make a global return to 
that region, tangencies between the different invariant manifolds play a significant 
role in the bifurcations of MMOs. We have not considered here secondary tangencies 
of the unstable manifold and the repelling slow manifold that are analogous to the 
bifurcations of secondary canards of a folded node analyzed by Wechselberger [32] ■ 

The second aspect of singular Hopf bifurcation that calls for further investigation 
is the presence of highly degenerate bifurcations within the normal form. In particular, 
we have seen that the system is very degenerate when 6 = 0. In analogy with the 
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study of codimension two bifurcations of equilibria |15j , the addition of higher order 
terms could resolve this degeneracy. We do not study that possibility here. Even in 
systems with a single time scale, the codimension three bifurcations that correspond 
to degenerate zero-Hopf bifurcations have not been studied thoroughly. 

Third, our investigation of the the points P where period-doubling bifurcation, 
fold of periodic orbits bifurcation and the slow flow has singular cycles meet at the 
degenerate homoclinic bifurcation remains incomplete. These points appear to be 
homoclinic orbit flips or inclination flips |2U|, but we have not reconciled aspects of 
the behavior found here with that predicted by the unfoldings of generic orbit flips 
or inclination flips. Investigations of phase portraits suggest that if B < 0, there is a 
homoclinic orbit flip at approximately {fi,A) — (i?/4, (— C± -y/C^ — 4i?)/2), provided 
that the square root evaluates to a real number. However, we find that the relative 
positions of the bifurcation curves are interchanged from those predicted by the the- 
ory, cf. [2 5) and figure 6.1 If i? > 0, we surmise the existence of an homoclinic orbit 




-3.5 -3 

H + 0.00248 



Fig. 6.1. Bifurcation diagram showing bifurcation curves of periodic orbits in the vicinity of 
homoclinic orbit flip at {fj,,A,B,C) = (-0.00248364,0.161803,-0.01,-0.1). Numbers with super- 
scripts s or u indicate the number of loops as well as the stability of periodic orbits in each region 
of parameter space. Unstable periodic orbits are of saddle-type. The labels H^ and H^ distinguish 
homoclinic bifurcation loci of periodic orbits with one and two loops respectively. 

[inclination] flip of type C [25] at approximately {p, A) = {B/A, {-C -VC^ - 4S)/2) 
if the parameter C is positive [negative], and at {fi, A) = (-B/4, (— C -|- \/C^ — 4S)/2) 
if parameter C is negative [positive], cf. section 3. 
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Appendix A. Bifurcation labels. 



Bifurcation label 


Bifurcation / explanation 


Codimension 


J^sup 


supercritical Hopf bifurcation 




Hsub 


subcritical Hopf bifurcation 




SN 


saddle-node bifurcation 




PD 


period-doubling bifurcation 




NS 


torus bifurcation / Neimark-Sacker bifurcation 




LPC 


fold of periodic orbits bifurcation / limit point of cycles bifurcation 




H 


homoclinic bifurcation 




T 


tangency of invariant manifolds bifurcation 




S 


approximate symmetry of the slow flow 




ZH 


zero Hopf bifurcation / fold Hopf bifurcation 


2 


GH 


generalized Hopf bifurcation / Bautin bifurcation / Gavrilov-Guckenheimer bifurcation 


2 


P 


see sections 4 and 5 


2 


Rl 


1:1 resonance on a torus bifurcation curve 


2 


R2 


1:2 resonance on a torus bifurcation curve 


2 


R3 


1:1 resonance on a torus bifurcation curve 


2 


R4 


1:4 resonance on a torus bifurcation curve 


2 


E 


fold of generalized Hopf bifurcations 


3 


F 


endpoint of the tangency curve crosses the torus bifurcation curve 


3 


G 


tangency curve coalesces with the Hopf bifurcation curve 


3 


H 


torus bifurcation curve intersect the Hopf bifurcation curve tangentially 


3 


I 


period-doubling bifurcation curve intersects the Hopf bifurcation curve tangentially 


3 


J 


fold of the curves where the slow-flow has singular cycles 


3 



Table A.l 

Table showing the abbreviations for bifurcations used in this paper, together with the full name or 
description of the bifurcation, as well as their codimension. Different names for the same bifurcation 
are separated by slashes. 
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Appendix B. Two dimensional bifurcation diagrams. 

This appendix is a catalog of {fi, A) bifurcation diagrams for system 



(1.2), 



eluding at least one sample diagram for each of the (B, C) regions la to Villa. 
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i^-TlPC^ ^ 
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Fig. B.l. Region la: {fJ., A) bifurcation diagram for {B,C) = (—0.01,0.1). This diagram shows 
homoclinic bifurcations as cyan curves. The top and bottom homoclinic curves have equilibria close 
the origin, while the middle two have equilibria distant from the origin. The first and the third 
homoclinic curves connect for 0{e^'^) values of ^, as do the second and fourth homoclinic curves. 
We conjecture that the endpoint of the top curve is close to where the periodic orbit born in the 
singular Hopf bifurcation ceases to be the only uj limit set for trajectories in the two-dimensional 
unstable manifold of the singular Hopf equilibrium, and that the endpoint of the bottom curve is 
close to where the periodic orbit bom in the singular Hopf bifurcation ceases to be the only a limit 
set for trajectories in the two-dimensional stable manifold of the singular Hopf equilibrium. The 
location of a tangential intersection of Sa with W"{Ef) is drawn with a dashed red line. 
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Fig. B.2. Region la: {fJ., A) bifurcation diagram for (B,C) = (—0.01,0.25). Observe that the 
tangency curve in some regions nearly coincides with the period- doubling and the torus bifurcation 
curves. This is not the case for all {B, C) in region la. For B < with large absolute value or C 
less positive, the torus bifurcation curve coincides with with the tangency curve on shorter segments, 
or not at all. For B < with small absolute value, the tangency curve tends to separate from the 
period- doubling curve closer to the point labeled P. 



0.02 



-0.02 



< -0.08 
-0.1 
-0.12 
-0.14 
-0.16 
-0.18 



^^^^.^^^subcrit.Hopf 




















- 


LPC 


GH 


^ 


^ 


-...^ (+T) 








S 




- 








supercrit\ 
Hopf 

7 


) 


^ 


> 


R2 




NS 


_ 




GH 




7 


■^y-^i*!) 


^ 










LPC^ 


R1 
.^ \ 


~^a^' 


^ 


^^tl- 


--^^^ 












supercrit. ^^ 
Hopt^^-^^ 


--^>, 


^ 


-subefit 


-Hqp^ 


_ _ _ _ SN_ 










S 





-"?R= 








SN 













-2 



7 
x10"' 



Fig. B.3. Region Ila: {fi,A) bifurcation diagram for (B,C) = (0.001,0.16). Note that the 
tangency curve partly coincides with the period-doubling curve, cf figure \3.1\ This is typically 
observed for very large values of C , cf. figure \C.1\ 
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Fig. B.4. Region Ila: {fJ., A) bifurcation diagram for (B,C) = (0.001,0.0896). The generalized 
Hopf bifurcations are close together, and vanish at {B,C) = (0.001,0.0894). 




Fig. B.5. Region Ilia: (m. A) bifurcation diagram for {B,C) = (0.001,0.08). 
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Fig. B.6. Region IVa: (/i,A) bifurcation diagram for {B,C) = (0.001,0.077). The tangency 
curve and the torus bifurcation curve do not meet in this diagram. 




Fig. B.7. Region Va: (^, A) bifurcation diagram for {B,C) = (0.001,0.075). 
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subcrit. Hopf 




Fig. B.8. Region Va/VIa: {fJ,, A) bifurcation diagram for {B,C) = (0.001,0.074). Relative to 
figure \B.6\ the resonances have crossed the Hopf curve, and a second torus bifurcation curve with 
R2, R3 and R4 resonances appeared. 
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Fig. B.9. Region Via: (fi,A) bifurcation diagram for (B,C) = (0.001,0.073). Although the 
torus bifurcation curve and period- doubling curve are at times very close, we did not detect inter- 
sections of the two curves. 
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supercrit. Hopf 



Fig. B.IO. Region Vila: {fi, A) bifurcation diagram for (B, C) = (0.001, 0.065). In this diagram, 
the curve LrgQ\, explained in section 4 and not included in other {lJ,,A) bifurcation diagrams, is 
plotted with a dashed black line. The inset shows that the curve lies on the branch of the periodic 
orbit that has real multipliers. Note that while the singular Hopf periodic orbit does not undergo 
any torus bifurcations for these values of B and C , there may be other {p., A) bifurcation diagrams 
in region IVa with torus bifurcations and the same kinds of codimension 2 bifurcations of periodic 
orbits that occur in region Ilia. 
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Fig. B.ll. Region Villa: (/i, A) bifurcation diagram for (B,C) = (0.02,0.1). 
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Appendix C. Locations of non-transverse two parameter dimensional 

slices. 



o T detaches from PD curve 
— G — GH bifurcations merge (E) 

o T detaches from NS curve (F) 
■ ■ o T attaches to Hopf curve (G) 

- e - R2 crosses Hopf curve 

- e - NS crosses Hopf curve (H) 

- e - PD crosses Hopf curve (i) 
— e — Singuiar cycies merge (J) 



CI -0.05 



e o . 



-2.5 -2 

log10(B) 



Fig. C.l. Approximate positions of degeneracies of {fi,A) bifurcation diagrams in (B,C) pa- 
rameter space for B,C' > 0. The legend lists eight types of degeneracies where {fJ.,A) slices of the 
parameter space fail to be transverse to codimension two bifurcation surfaces. 



Inspecting {fJ-, A) bifurcation diagrams for a large number of parameters {B,C), 
it is easy to see that the catalog of diagrams completely describes the bifurcations 
of the types considered in this paper for B < 0. The completeness of the diagrams 
for B > was verified by computing curves that separate the {B, C) plane into re- 



gions with qualitatively different (/i, A) bifurcation diagrams as shown in figure C.l 



Where applicable, the curve labels introduced in section 4 are included in parentheses. 
For numerical reasons, \B\ was restricted to lie between 10^** and 10^'^. The positions 
of the folds of the generalized Hopf bifurcations and of the points where the slow flow 
has singular cycles are calculated using asymptotic relations. All other data points 
were obtained from sequences of (/i. A) bifurcation diagrams. The error in the ordi- 
nates of the data points is significant, but the order of the data points on the ordinate 
is correct for each value of B. Note how the dashed green line intersects several other 
curves, corresponding to further degeneracy in the set of (/i. A) bifurcation diagrams. 



Figure B.2 illustrates the degeneracies that can occur in (/i,v4) bifurcation di- 
agrams if C is too large relative to jSp^. The saddle-node bifurcation is close to 
intersecting the dashed curve where the slow-flow has singular cycles. 
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Appendix D. Numerical methods for the computation of the tangency 
curve. 

This appendix describes the two mimerical methods used for the computation of 
the tangency curve. The first method is a shooting method, used by Desroches et al 
([B]) to compute the tangency curve for the Koper model discussed in section 5. The 
second method uses collocation and the software package AUTO [H]. 

The tangency bifurcation marks the onset of a bistability of trajectories in the 
unstable manifold W'^{Ef) of the saddle- focus Ef near the origin: the repelling slow 
manifold Sr separates trajectories that flow to A = — cx) from those that remain in the 
fold region. This observation motivates the use of a shooting algorithm to compute the 
position of the tangency curve in parameter space: given a set of parameters, a grid 
of initial conditions in a linear approximation of W"{Ef) is integrated numerically for 
a long time interval. The tangency curve separates the parameters where at least one 
trajectory reaches A = — oo from the parameters where all trajectories approach a 
bounded attractor. This curve can be computed by a 2-parameter predictor-corrector 
continuation method starting from an initial parameter on the tangency curve, where 
the correction may be implemented using interval bisection. 

The main disadvantage of this method is that it requires choices regarding the 
number of grid points and, more importantly, the length of the integration-time. 
Unless convergence to an attractor such as a limit cycle can be determined easily, 
this method cannot determine in an automated fashion whether a trajectory that 
has stayed in the fold region for a long period of time will remain in the fold region 
forever or leave it after a phase of transient behavior. For example, if Ef is close to 
a Hopf bifurcation curve, and the real part of its eigenvalues are thus close to zero, 
convergence to the periodic orbit is slow and takes a very long time. In this case, the 
algorithm may give incorrect results unless the integration time is chosen adaptively 

The run-time of the algorithm described above is proportional to the number of 
grid points and the integration-time, as the bulk of the computation time of the al- 
gorithm is spent on numerical time-integration of trajectories that remain in the fold 
region for the entire integration time. Slight speed-ups can be achieved by choosing 
the order of time integrations for grids of initial conditions to start with points that 
are far apart in the fundamental domain. Sometimes, this locates a trajectory that 
escapes the fold region more quickly than choosing initial points consecutively along 
a segment of a fundamental domain of the linearized unstable manifold. We found 
that using more than 10 grid points did not improve the accuracy of the method 
significantly. Larger speed-ups are possible by parallelizing the algorithm. Since the 
tangency curve is smooth, parallelism can be implemented using a two-pass approach: 
in the first pass, a coarse approximation of the position of the curve is computed on 
a single processor with the algorithm described above, using large step-sizes, few grid 
points, and thus little computation time. The second pass interpolates between data 
points of the first pass, using larger numbers of interpolated points, more grid points 
and long integration times on multiple processors. This second pass is embarrassingly 
parallel. 

The tangency curve can also be computed via continuation of a boundary value 
problem in AUTO. The boundary conditions require that trajectories begin on a fixed 
ray in the linear approximation of the equilibrium's unstable manifold and end on the 
parabola Y = A^ -I- 5. An initial trajectory segment satisfying these requirements, 
obtained for example using the first method, can then be continued in one parameter. 
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AUTO's fold detectfon determines when the continued solution folds, i.e. W"{Ef) 
and Sr intersect tangentially. Switching to a continuation of the fold, the tangency 
can now be continued in two parameters. Computing the tangency curve in AUTO is 
in general much faster than the shooting method, but the shooting method provides 
additional insight into the fate of the trajectories in W'^{Ef) before and after the 
tangency. As with the shooting method, the boundary value algorithm breaks down 
close to the Hopf curve, unless the number of mesh points is increased appropriately. 
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